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Heat-Capacity Measurements of Energy-Gap Nodes of the Heavy-Fermion 
Superconductor Celrlns Deep inside the Pressure-Dependent Dome Structure of its 

Superconducting Phase Diagram 
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We use heat capacity measurements as a function of field rotation to identify the nodal gap 
structure of Celrlns at pressures to 2.05 GPa, deep inside its superconducting dome. A four-fold 
oscillation in the heat capacity at 0.3 K is observed for all pressures but with its sign reversed 
between 1.50 and 0.90 GPa. On the basis of recent theoretical models for the field-angle dependent 
specific heat, all data, including the sign reversal, imply a d^2_y2 order parameter with nodes along 
[110], which constrains theoretical models of the pairing mechanism in Celrlns. 

PACS numbers: 74.20. Rp, 74.25.Bt, 74.40. Kb, 74.70.Tx 



The Ce-based heavy- fermion family CeMIn5(M=Rh, 
Co, Ir) has attracted attention due to their exotic prop- 
erties and interesting interplay between magnetism and 
superconductivity CcRhlns is antiferromagnctic 

(AFM) at ambient pressure but becomes superconduct- 
ing as pressure suppresses the AFM order and tunes 
CeRhlns to a quantum-critical point (QCP) 0, indicat- 
ing an intimate relationship between AFM fluctuations 
and unconventional superconductivity. CeCoIns, with 
the highest superconducting Tc=2.3 K among Cc-bascd 
heavy- fermion materials 0], is a nodal superconductor 
with a dx2_y2 order parameter (OP) that is revealed con- 
sistently in a variety of exprimental measurements 043 ■ 
Cd or Hg doping, acting as an effective negative pressure, 
tunes the ground state of CeCoIns from superconducting 
to AFM [a [13. The resulting temperature-doping phase 
diagram is very similar to the temperature-pressure dia- 
gram of CeRhlns , revealing a common response of these 
materials to the interplay between magnetism and super- 
conductivity. Celrlns is also superconducting at ambient 
pressure with bulk Tc ^ 0.4 K but resistive Tc ^ 1.2 K 
[lll . [T^ . This difference in TcS is found as well when mag- 
netic order coexists with superconductivity in CeRhlns 
and doped CeCoIns [l^ , but there in no long-range AFM 
in Celrlns. Further, the response of Celrlns to effective 
pressure is different, as shown in the inset of Figllja). 
Rh-doping in CeRh^^Iri-^jIns acts as an effective nega- 
tive pressure and induces a dome of superconductivity in 
proximity to AFM at x>0.5. Tc goes to zero or approx- 
imately zero at x~ 0.1 before it rises again with applied 
pressure to form a second dome [13, ■ This diagram is 
very similar to that of CeCu2(Sii_a;Gea;)2 as a function of 
pressure [3 and has led to the suggestion that the sec- 
ond superconducting dome (SCII), so distant from AFM 
order and thus less associated with AFM flutuations 17| , 



may be due to another pairing mechanism, such as Ce- 
valence fluctuations 18[, and possibly supports a differ- 
ent OP symmetry. 

Identifying the nodal gap structure and pairing sym- 



metry in Celrlns is an important step toward resolv- 
ing the pairing mechanism. Power laws in heat capac- 
ity, thermal conductivity [l^ and penetration depth 
arc consistent with line nodes in the gap, but these do 
not probe the node positions, which is needed to de- 
termine the allowed OP symmetry. A large anisotropy 
between in-plane and out-of-plane thermal conductivity 
measurements on Celrlns has led to a proposed hybrid 
gap, kz{kx+iky) [21, 22|. In contrast, thermal conductiv- 
ity studies in a magnetic fleld that is rotated within the 
basal plane imply a (1^2 _y2 gap with line nodes along the 
[110] direction [isl. similar to CeCoIns [^ and CeRhlns 
under pressure [24| . This discrepancy may be due in 
part to both measurements being made at ambient pres- 
sure where the superconductivity of Celrlns sits right at 
the margin of the two superconducting domes and may 
be influenced by the competition between residual AFM 
fluctuations and other pairing interactions, if any. Prob- 
ing the gap structure in the superconducting state deep 
inside SCII avoids these complications. 

In this Letter, we report heat capacity measurements 
of Celrlns under pressures up to 2.05 GPa as a mag- 
netic fleld is rotated around and through its tetragonal 
c-axis. A four-fold oscillation of the in-plane heat capac- 
ity indicates a gap of dx2-y2 symmetry with nodes along 
[110], even deep inside SCII. This OP symmetry should 
be captured in mechanistic models of superconductivity 
of Celrlns in SCII. 

Celrlns single crystals were grown out of In-rich flux 
and screened by magnetic susceptibility to ensure the 
absence of free In. The sample was mounted in a Be- 
Cu/NiCrAl hybrid clamp- type pressure cell with silicone 
fluid as the pressure-transmitting medium, which pro- 
vides a very nearly hydrostatic environment to 25 kbar. 
The resistive superconducting transition temperature of 
Sn was measured to determine the pressure at low tem- 
peratures. The heat capacity of Celrlns Cac under pres- 
sure was measured by an ac calorimetric method wherein 
heating from an ac current generates an oscillation in the 
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sample temperature AT^c that is measured by a field- 
calibrated Cr/AuFe(0.07 %) thermocouple through its ac 
response; the heat capacity is inversely proportional to 
ATac, i.e., Cac oc 1/ATac- Magnetic field rotation was 
provided by a triple-axis vector magnet that accomodates 
a ^He cryostat with the pressure cell inside. 
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FIG. 1: (color online), (a) Temperature-dependent heat ca- 
pacity of Celrlns at pressures P=2.05 (circles), 1.50 (dia- 
monds) and 0.90 GPa (stars) normalized to the values at 
T=1.8 K. The inset is the phase diagram of Ce(Ir,Rh)In5 as 
a function of Rh doping and pressure, taken from Ref[2^; (b) 
Heat capacity of Celrlns at T=0.3 K and P=2.05 GPa as a 
function of magnetic field applied along [001] (squares), [100] 
(diamonds) and [110] (circles) direction, respectively. 

Figure [I{a) shows the temperature dependence of the 
heat capacity Cac under pressures P=0.9, 1.5 and 2.05 
GPa going deeper into SCII with increased bulk Tc=0.50, 
0.75 and 0.85 K, respectively. The magnetic field de- 
pendent heat capacity Cac is shown in Fig. [Ijb) for 
Celrlns at 0.3 K and 2.05 GPa with the field applied 
along three major crystallographic orientations: [001], 
[100] and [110]. The upper critical field along the c-axis 
H^2 is smaller than that within the ab plane by a factor of 



2, similar to the case of CeCoIus [26[ and CeRhlns [24 1. 
Within the ab plane, h[.2°''' along [100] is reproducbily a 
little larger than hJ.j^''^ in the [110] direction, probably 
refiecting the d-wave superconducting OP symmetry as 
discussed later. 

Figure [2] shows the field-angle heat capacity of Celrln.g 
at 2.05 GPa deep inside SCII with magnetic field rotated 
in the ab (ac) plane. For the field rotated within the ab 
plane, Cac/T is featureless in the normal state at 1.8 K 
even with magnetic field up to 8.0 kG. An apparent four- 



fold oscillation is present in ^fr- (0) when cooling down to 
0.3 K in the superconducting state and applying field at 
2.0 kG. A polar sweep at 0.3 K with H=2.0 kG is shown 
in Fig. [2]^b) and a two-fold modulation is observed, which 
can be fitted to the function Cac[0) = Co+C2\cos9\, most 
likely caused by the Hc2 anisotropy between c-axis and 
ab plane. 
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FIG. 2: (color online). Angle-dependent heat capacity of 
Celrlns at 2.05 GPa. (a) Azimuthal rotation (<^) of the 
magnetic field in the ab-plane in the normal (T=1.8 K and 
H=8.0 kG) and superconducting state (T=0.3 K and H=2.0 
kG). The solid line is a fit to a four- fold modulation model 
Cac = Co + C2Coa2{cj> — 4>o) -f CiCosA4>. (b) Polar rotation 
[6) of the magnetic field (2.0 kG) in the ac-plane of Celrlns 
at T=0.3 K. The solid line is a fit to a two-fold oscillation 
model Cac=Co-|-C2|cos6|, where the modulation arises from 
IIc2 anisotropy. 

Theoretical and experimental progress has established 
field-angle dependent heat capacity (and thermal con- 
ductivity) as a powerful probe of nodal locations in the 
superconducting gap structure d, [13, [13, HI]. In the 
Abrikosov state where supercurrents circulate around the 
vortex core, extended quasiparticles (QPs) experience a 
Doppler shift, 5E ~ pp • Vg, where pf is the Fermi mo- 
mentum of QPs and Vg is the supercurrent velocity. In 
nodal superconductors at low temperatures, the Doppler 
energy shift plays an important role in the QP excita- 
tions. Near the nodes, the local superconducting gap is 
sufficiently small that the Doppler shift breaks Cooper 
pairs, A(pf) < |pf ■ Vs|, and, consequently, the QP den- 
sity of states (DOS) depends on the field orientation rel- 
ative to the position of gap nodes. In a classical picture, 
the QP DOS and thus heat capacity is a minimum for the 
field along nodal directions and a maximum along antin- 
odal directions. In the case of a d-wave superconductor, a 
four-fold oscillation emerges in the superconducting state 
when the field is rotated in the basal plane and the local 
heat capacity minimum refiects a nodal position on the 
Fermi surface [i^. The four-fold oscillation observed in 
an azimuthal field sweep for Celrlns at 2.05 GPa implies 
the presence of vertical line nodes perpendicular to the 
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ab-plane as expected for a gap with d-wave symmetry. A 
hybrid gap with horizontal hne nodes in the equator is 
inconsistent with the four-fold oscillation observed here. 
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FIG. 3: (color onlme). (a)-(d) Heat capacity of Celrliis at 
T=0.3 K and P=2.05 GPa as a function of azimuthal field 
rotation for H=2.0, 5.0, 8.0 and 9.0 kG, respectively. The 
upper critical field along [100] is f 1.0 kG at this temperature 
and pressure, (e-g) Azimuthal field dependence of the heat 
capacity of Celrlus at T=0.3 K and P=0.90 GPa for H=2.0, 
3.5 and 5.0 kG, respectively. The upper critical field along 
[100] is 6.70 kG. Solid lines are a fit to the four-fold oscil- 
lation model and the vertical dashed line indicates the [110] 
direction. The two-fold and four-fold components of the fit 
in (a)-(c) and (e)-(f) are plotted as dotted and dashed lines, 
respectively, (h) The four-fold oscillation component C40 at 
T=0.3 K as a function of H/hI.j'"'' normalized to the corre- 
sponding ACaciH) = C{H)-C{0) at P=2.05 (squares), 1.50 
(diamonds) and 0.90 GPa (circles), respectively. 

The angular location of maxima and minima in Cac((0) 
allows a distinction in the OP symmetry between dxy and 
d^2_y2. Figs. El (a)-(d) show the heat capacity Cac of 
Cclrlns at 2.05 GPa and 0.3 K as a function of field angle 
(j> against the [100] direction in the ab plane at different 
fields up to 9.5 kG. Data at 9.5 kG (Fig. |3](d)) are qual- 
itatively different from those at lower fields (Figs. [3][a)- 
(c)): at 9.5 kG local minima are shifted by 45 degrees and 
the oscillation amplitude is one order of magnitude larger 
relative to those in (a)-(c). We note that this magnetic 
field is very close to Kc2 for [100] and [110] directions. 
The in-plane Hc2 anisotropy will contribute four-fold os- 
cillation in Cac(0) with local minina in [110] directions, 
providing an indirect way to check the crystalline orienta- 
tion. The offset of C^^ between [110] and [110] in FigEJd) 
indicates a small sample misalignment between the crys- 
tal ab-planc and the field xy-planc, which introduces a 



two-fold oscillation in the field rotation heat capacity. 
The total heat capacity in magnetic field can thus be 
writtern as Cac(0)=Co+C2Cos2((/) — (/)o)-l-C4COs4(/), where 
Co is the zero-field heat capacity coming from thermally 
excited QPs and phonons, the two- fold term is due to 
the sample misalignment, and the four- fold oscillation 
C4Cos4(/) arises from nodal structure of the superconduct- 
ing OP. 

Though the local minima of Cac{(t>) are all in the [100] 
direction and maxima are along [110] for different mag- 
netic fields as shown in Fig. |3Ia)-(c), the naive classical 
picture would imply a d^y OP symmetry with line nodes 
along [100] directions for Celrlus at 2.05 GPa. However, 
Vorontsov and Vekhter have argued that the classical pic- 
ture holds only when T/Tc <C 1 and H/Hc2 <C 1, and, 
consequently, the angular dependence of four-fold oscil- 
lations shows a complex evolution across the whole H-T 
phase diagram [s^ . In general, anisotropic scattering of 
QPs due to vortices produced by the magnetic field plays 
an important role in the DOS variation and thus the heat 
capacity as a function of the field orientation, which is 
missed in the classical picture. At an intermediate tem- 
perature or moderately high field, angular anisotropy in 
the heat capacity changes sign and the four-fold oscil- 
lation is inverted, showing maxima rather than minima 
when the field is along the nodal direction. According 
to their calculation of the field-angle dependent heat ca- 
pacity of CeCoIns, where a corrugated cylindrical FS is 
assumed, our field rotation heat capacity measurments at 
0.3 K, with T/Tc ^ 0.35, are in the intermediate region of 
the H-T phase diagram and a four-fold component with 
maxima in the nodal directions is expected. Our observa- 
tion of oscillation maxima in [110] directions for CcIrlus 
at 2.05 GPa indicates that the line nodes are along di- 
agonals of the crystal lattice and favors dx2-y2 symme- 
try over dxy for Celrlus at 2.05GPa deep inside SCII. 
Similar behaviors observed for Celrlus at 1.50 GPa with 
T/Tc ^0.40 (not shown here) also are consistent with a 
dx2-y2 OP symmetry. 

This interpretation is supported by results plotted in 
Figs. [3] (e)-(g) that show the field-rotation heat capac- 
ity Cac of Celrlns at 0.90 GPa and 0.3 K as a function 
of field angle (j) for H=2.0, 3.5 and 5.0 kC, respectively. 
In contrast to the data at 2.05 GPa, four-fold oscilla- 
tions at 2.0 and 3.5 kG are inverted and show minima 
along [110]. Figure |3] (h) is a plot of the normalized 
four- fold oscillation amplitude C4/AC(H) as a function 
of normalized field H/Hc at 0.3 K for Celrlus at 2.05, 1.50 
and 0.90 GPa, where AC(H)=C(H)-C(0) is the change of 
heat capacity in magnetic field. With similar amplitudes 
(3 - 6%), the four- fold oscillations have the same sign 
at 2.05 and 1.50 GPa but a different sign at 0.90 GPa. 
We summarize all of our results at different pressures 
and magnetic fields and compare them to the theoretical 
expectations (3l| in Fig. |4] . This comparison assumes 
the similarity between CeCoIns and Celrlns mentioned 
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FIG. 4: (color online). The calculated four- fold oscillation of 
field-rotation heat capacity in a normalized H-T phase dia- 
gram, taken from Ref. l3ll . Our data points are plotted in 
the H-T phase diagram for Celrlns at 2.05 (diamonds), 1.50 
(stars) and 0.90 GPa (circles). 

above. While the measurements at 2.05 and 1.50 GPa 
are located in the intermediate region, the points at 0.90 
GPa and T/Tc ~ 0.6 are near the boundary of the low 
H/Hc, high T/Tc shaded region where the four-fold os- 
cillation swtiches its sign again from that expected in the 
intermediate region and exhibit minima for field in the 
nodal directions. We note that the detailed division of 
different regions depends on the exact shape of the Fermi 
surface, and the general trend of sign change from inter- 
mediate to high temperature ranges should survive [sij. 
Consequently, oscillation minima along [110] for Celrlns 
at 0.90 GPa also are consistent with dx2-y2 OP symme- 
try. 

An interesting point is that we observe an in-plane up- 
per critical field anisotropy, with Hc2 in the [110] direc- 
tion always smaller than along [100] for all p ressures.On 
the basis of model calculations in Ref. ISOj, this upper 
critical field anisotropy between nodal and antinodal di- 
rection appears naturally as a result of the d-wave gap 
symmetry with W^f'^ ^^anUnode ^^so stress that 

the in-plane Hc2 anisotropy cannot be the origin of the 
observed four-fold heat capacity oscillation because the 
oscillation change its sign between 2.05 and 0.90 GPa but 
the in-plane Hc2 anisotropy remains the same. Wc thus 
conclude that the line nodes are in [110] directions and 
the gap symmetry of Celrlns stays dx-2-y2 across SCII. 

Because of proximity to AFM order, it is natural to 
consider that the OP symmetry in SCI is Ax-2_y-2 as it is in 
CeCoIns or CeRhln.g under pressure, which also are near 
AFM order. On the other hand, our field-rotation heat 
capacity results also favor dx^-yi. gap symmetry deep in- 
side SCII, up to 2.05 GPa. In this regard, recent trans- 
port measurements observe non-Fermi liquid behaviors 
and claim that AFM fluctuations survive deep into SCII 
psj . It is then an interesting question why the supercon- 
ducting phase apparently is divided into two regions but 



with the same gap symmetry. One possible scenario is 
that there is some hidden competing order that destroys 
bulk superconductivity in Celrlns near ambient pressure 
and produces a minimum in the bulk superconducting 
transition temperature. More studies are required to 
clarify the origin of two superconducting domes in the 
Celrlns system. 
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